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Abstract 

Using the heat kernel method and the analytic continuation of the zeta function, we calculate 
the canonical and improved vacuum stress tensors, {T^j_y{x)) and {Q^p{x)), associated with a 
massless scalar field confined in the interior of an infinitely long rectangular waveguide. The 
local depence of the renormalized energy for two special configurations when the total energy is 
positive and negative are presented using {Tqq{x)) and (6oo(^))- From the stress tensors we obtain 
the local Casimir forces in all walls by introducing a particular external configuration. It is shown 
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that this external configuration cannot give account of the edge divergences of the local forces. 
The local form of the forces is obtained for three special configurations. 



1 



1 Introduction 



In a previous paper, the one-loop renormalization of the anisotropic scalar model was performed, 
assuming that the fields were defined in a d-dimensional Euclidean space where the first d — 1 
coordinates are unbounded, while the last one lies in the interval [0,L] |T|]. The authors analysed 
the vacuum activity of massive scalar fields assuming different boundary conditions on the plates, 
namely Dirichlet-Dirichlet {D — D) boundary conditions and also Neumann- Neumann (A^ — A^) 
boundary conditions. They obtained two different results. The first one has been obtained pre- 
viously by many authors, and is the fact that to renormalize the theory we have to introduce 
counterterms as surface interactions. The second one is the fact that the tadpole graph for DD 
and for A^A^ have the same z dependent part in modulus but with opposite signs. This second 
result has been obtained by DeWitt and also Deutsch and Candelas 0. In ref.[|l|, the authors 
also investigated the relevance of this fact to eliminate the surface divergences. 

More than thirty years ago, the local version of the Casimir original problem was performed 
by Brown and Maclay They obtained a constant stress-energy tensor due to the cancelation 
of the electric and magnetic sectors, showing the uniformity of the vacuum of the electromagnetic 
field for this configuration. The result is due in part to the particular field involved and also 
to the simplicity of the parallel plane geometry. A similar cancelation can be arranged for the 
scalar field by computing the improved stress-energy tensor, but in a more complicated rectangular 
geometry with the presence of edges and corners (e.g. within a rectangular waveguide), we expect 
an answer strongly non-uniform. The aim of this paper is to generalize part of the results of 
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Fosco and Svaiter introducing edges in the domain where the fields are defined, calculating 
the renormalized stress-energy tensor of a massless scalar field in an infinitely long rectangular 
waveguide. As was stressed by Maclay, geometries with corners present special problems with 
respect to vacuum energy; nevertheless these issues have received little attention in the literature 
P]. It is commonly accepted that, the understanding of the renormalization of the stress-energy 
tensor of quantum fields in the presence of classical boundaries should throw light on the more 
difficult case, where there is the added complication of local curvature effects [Q. 

It is well known that there are two quantities which might be expected to correspond to the 
total renormalized energy of quantum fields 0. The first one is called the mode sum energy and 
its definition is 

I mode /" • 1 



{EZT = l du;-uj[N{u;)-No{u;)], (1) 

where is the zero point energy for each mode, N{uj)duj is the number of modes with frequencies 
between u and u + duj m. the presence of boundaries and NQ{ijj)duj is the corresponding quantity 
evaluated in empty space. The above equation gives the renormalized sum of the zero point energy 
for each mode. The second one is the volume integral of the renormalized energy density (-E)™^ 
obtained by the Green's function method [§,3- A recent investigation of {E)^"^ in rectangular 



geometries was given by Svaiter and colaborators |jTO|, [III. A seminal paper studying this kind 



of geometry was made by Ambjorn and Wolfram [12], and more recently Milton and Ng studied 



the Casimir effect in (2 + 1) Maxwell-Chern-Simons electrodynamics in a rectangular domain 



131 . Since these definitions deal with integrated quantities, surface divergence problems do not 



appear in the calculations. Although global effects are more accessible to experiments, it is quite 
important to understand how the global effect is obtained from the local version. This issue has 
recently been studied by Actor and Bender 

In ref. ^ the author studied the use of the zeta function method to find the effective action 
associated with a scalar field defined in the interior of the infinitely long waveguide, while in ref. 
the authors use the same method to compute the stress-energy tensor for various rectangular 
geometries. Using the relation between the local force density and the discontinuity of the stress- 
energy tensor across the boundaries, they computed the local Casimir forces, which exhibited 
strong position dependence. 

In this paper we are interested in calculating local quantities in the presence of surfaces and 



edges. As was stressed by Dowker and Kennedy [14| and also Actor and Bender M, to study the 



local problem in the infinitely long rectangular waveguide, it is necessary to present the local form 
of the analytic continuation of the local zeta function in the rectangle. Note that our choice of a 
rectangular cavity is related to the fact that the modes of the field in this geometric configuration 
are well known and an exact calculation can be done. 

The organization of the paper is the following: In section II a brief review of the zeta function 
method is presented. In section III we calculate the vacuum expectation value of the canonical 
and improved stress-energy tensors associated with a massless scalar field using the zeta function 
method in the infinitely long rectangular waveguide. We also show in section III the relation that 
exists between the local version and the global version of the Casimir energy for the waveguide. 
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In section IV, we use the results of the previous section to compute the local forces. In order to 
do this we introduce an external configuration (such that the interior region is the waveguide) 
for which the components of the stress-tensors are known everywhere. Conclusions are given in 
section V. In this paper we use h — c— 1. 



2 The canonical and improved stress tensors and the zeta 
function method 

In this section we will describe the basic procedure to compute the renormalized vacuum expec- 
tation value of the stress-energy tensor for a real scalar field. Our aproach will be based on the 
zeta function method. 

For a real scalar field defined in a four dimensional spacetime, distorted by static boundaries, 
we can use the Fourier standard expansion 

<P{x) = Yl ke--°-"0„(f) + at e-°-"0:(f)] . (2) 

Assuming that the manifold is static, i.e., that it possesses a timelike Killing vector field, it is 
possible to show that there is a complete set of spatial modes {(t)n{x)} satisfying a Schrodinger- 
like equation 

- A(f)n{x) =Ujl(l)n{x), (3) 

where over these modes we will impose certain boundary conditions. Here we are concerned only 
with Dirichlet boundary conditions, although the generalization to Neumann boundary conditions 



is straightfoward. Since the set of modes are orthonormal and complete, one then readily 

verifies that the equal-time canonical commutation relations imply the usual commutation relation 
between annihilation and creation operators of the quanta of the field. 

The main point of interest for us will be the renormalized stress-energy tensor of the scalar field 
confined in the interior of the rectangular infinitely long waveguide. The canonical and improved 
stress tensors of a real massless scalar field are given by 



and 



T^,u{x) = ^ [d^(f> dy(j) + ducj) df,(f) - r]f,yda 



(4) 



(5) 



where r/^jy is the Minkowski metric. One way to write the vacuum expectation value of T^y(x) 
using eq.(0) is 



{T^,y{x)) = lim - 



d d d d 
+ 



d d 



(mm) 



(6) 



dx^^ dy^ dx^ dy^^ dx" dy, 

where {(j){x)(j){y)) is the vacuum expectation value of the product of the fields in two different 
points. (An equivalent relation exists for (0^,^).) Using the commutation relations between anni- 
hilation and creation operators, the quantity in eq.(|D, can be written as 



= XI TT- exp(-i (xo - yo) ^„)0„(f )0* 



(7) 



It is clear that (T^,y(x)) can be obtained from the bilocal sum given by eq.(|^). The bilocal (spectral) 
sum in eq.(|^) diverges and needs a regularization and renormalization procedure. A convenient 



method is to set xq = yo and replace uj~^ in eq.(0) by cj^^'' with s complex, initially holding for 
Re{s) > and sufficiently large to guarantee convergence even for x = y , followed by analytic 
continuation in s. 

Let us work with a compact manifold M with or without boundaries. The diagonal zeta 
function associated with some elliptic, semi-positive and self-adjoint differential operator D will 
be defined by ({s\D). Let 0„(a;) and be the spectral decomposition of i5 in a complete normal 
set of eigenf unctions 0„(x) with eigenvalues An, i.e. 

D(t)n{x) = Xn(pn{x) (8) 

where (f)n{x) = {x\n). Since the eigenfuctions (j)n{x) form a complete and normal set it is possible 
to define the generalized zeta operator associated with D as 

C(.|D)=/.-t^> (9) 

where we introduce the parameter /i with dimensions of mass in order to have a dimensionless 
quantity raised to a complex power and the prime sign indicates that the zero eigenvalue of D 
must be ommited. The generalized zeta function associated with the operator D is defined by 

C{s\D) = I?' j^d-i{x)lx\D^'\x), (10) 

where d''){x) is the measure on M. We have then to consider the bilocal zeta function 

as I = /i2^^(c.^^)->„(f)0;(^-'), (11) 

n 
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which has abscissa of convergence Re{s) = |. Since the modes 4>n{x) form an orthonormal set 
then the passage from the local to the more familiar global zeta function is straightforward for 
Re{s) > |. This can be done integrating the bilocal zeta function, i.e. 

as)=^^'^ ^d^ix)as\x,x)=^^'^J:Hy\ Reis)>l (12) 

n ^ 

A careful analysis of the analytic extension of the global zeta function associated with some 
differential operator defined in compact manifold with or without boundaries can be found in ref. 



rsf . Going back to the local case in the analytic extension of the local zeta function to the whole 
complex plane (to the region Re{s) < |), it will appear poles related with the geometry of the 
manifold. For sake of simplicity we will omit the fi factor in the following. 

The function given by eq.(|lT]) is related to the heat kernel by a Mellin transform 



1 

C(s\x,y) = —- dtf-'K{t\x,y) (13) 
1 [s) Jo 

where 

K{t\x,y) = Y.e-'^^n{xmy) (14) 

n 

is the heat kernel satisfying the same boundary conditions that we choose to the complete set of 
modes 4>n{x). It is possible to express the vacuum expectation value of the canonical stress tensor 
given by eq.(|^) in terms of the modes and also the frequencies Un- We have 

(Tooix)) = ]Y.^n\Kf + (15) 
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and 



{Tii{x)) = ^J2ujn\(pnf -]j2-r ^ + ^Hrr 1^^*^^'^ 



4 



2 1 

2 



i not summed. It is easy to see that for 0„ real and for plane waves 



(16) 



(Toiix)) = 



(17) 



and finally 



For the improved stress tensor we have: 



5 11^ 



(18) 



(19) 



(20) 

i not summed, 

(60.(0;)) = (To,(.t)) = 0, (21) 



1 E - ['/'n + ( Wn) ^T^J- (22) 



In the next section we will identify the divergences and the finite parts that appear in the 
vacuum expectation value of the canonical and the improved stress tensors of a real massless 
scalar field satisfying Dirichlet boundary conditions in all walls of an infinitely long rectangular 
waveguide. 
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3 Canonical and improved stress-energy tensor of a mass- 
less scalar field confined within a rectangular waveguide 

In this section we will apply the local zeta function method to calculate the renormalized vacuum 
expectation values of the canonical and improved stress-energy tensors of a massless scalar field 
confined within an infinitely long rectangular waveguide. Let the waveguide be oriented along the 
axis in such a way that the field is defined free in the region 

= X = (xi, X2, X3) : < Xi < a, < X2 < b C R^, (23) 

with Dirichlet boundary conditions at Xi = and xi = a and also X2 = and X2 = b. The spatial 
modes are given by: 

(X) = (-) s,u — s,u — e", (24) 



with mi,2 = 1, 2, 3, ... and —00 < < 00. The eigenvalues are given by 



where n denotes the collective indices {mi, 1712, k^)- Substituting eq.(p4D in eq.([T^) the heat-kernel 
can be written as: 



K{t\x,y) = ^e-*--(/)„(x)0™(y) 

/oo °° ( 

dks exp <^ -t 



1 / 4 \ roo 
27r \ab) 



mi,m2=l 



X sm( )sm( — ; )sm( ^) sm( — ^_)e««3ia;3 s/3j_ /26) 

a b a b 
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The free spacetime part can be integrated imediately: 



{x3 - ysY 
At 



(27) 



yielding 



K{t \ x,y) = — (47rt) 2 exp 
ab 



{x3 - ys 



i2n 



00 

X ^ exp 

mi=l 

00 

X ^ exp 

m2=l 



-t 



niiTT 



a 



/ 17127!- 



minxi miTryi, 
sml ) sml 



. ,m2'KX2. . ,m2Txy2. 
sm — : sm — : 



(28) 



Using trigonometric identities and also the Jacobi 6'-function identity 



exp(— m^x) cos{m2Tih) 



771=1 



{^n + hy — 

X 



(29) 



one finds to the heat-kernel 



K {t \ x,y) = (Airt) 2 exp 



2n 



00 t 
X ) exp 

ni=— 00 I 
00 I" 

X J2 ) exp 

n2=— 00 I 



(3^3 - 1/3) 

4t 



- [2nia + (xi - yi)] 

4t 

- [2n2b + {x2 - y2)Y 

At 



2n 



exp 



exp 



[2nia + (xi + yi)] 



2n 



it 



- [2n2b + {X2 + y2)] 
At 



2n 



■ (30) 



As we discussed before to find the bilocal zeta function we need to perform the Mellin transform 
of the heat-kernel given by eq.(|30|). All terms of eq.(^) can be integrated using [|l( 



5 A 3 

/ dtf-ieM ) = A^^^r(J-s). 

JQ t 2 



(31) 
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After a straightforward calculation we have 



c{s\x,y)= ^/;^3„;, E (^1 + ^2 + ^3 + ^4) 



(47r)2 r(s) ni,n2=-oo 

where = n2, a;, j = 1, 2, 3, 4, are given by 



nia + 



n,a+ ("^:^^^ l +f-2&+^^^^l +(^)^ 



(32) 



(33) 
(34) 
(35) 
(36) 



We see that divergences appear in the local zeta function C{s \ x,y) in the limit y ^ x. We note 



that ({s I X, x) has surface divergences when Re{s) < | . The term Z2 (0, 0, x, x) = {x2)'^'^~^, for 
example, diverges when X2 — in this case. 

In order to calculate the components of {T^j^y{x)) we have to evaluate the mode sums given by 



eqs.(|T5D-(p!8D. One then readily verifies that 



E |0n|^ = C(S = -- I f, X) 



;Fo(x) 



(37) 



where the expression for Fo{x) is given by 



Fo{S) = E liniay + in^b) 



-2 



+ 



ni,n2=— 00 



12 



(nifl)^ + (n2& + 
{riia + xi)^ + (^26)" 



-2 



+ 



+ 



+ 



The other terms that we need are given by: 



Substituting eg. (p2D in eq.(^) for i = 1, we have: 



2 d d 1 I ^ 

- hm — — C(s = - I x,y). 



1 



where the functions -Di(x) and Fi (x) are defined by 

o 

[(nia)^ + (n2& + ^2)^ [nia]^ + 
(raifl + xi) + (n26) [riia + xi] + 



ni ,n2=— 00 



+ 



(nia + xi) +(n2& + X2) [nio + Xi] 



-3 



and 



Fi{x) 



-2 



ni ,712=— 00 



- [(riitt)^ + {n2h + X2)' 
+ (riia + Xi)^ + (71.26)^ 



+ 



-2 



+ 



(nia + xif + (?22& + 3:2)' 



-2 
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J2 — 1^20^ 



47r2 



where the functions D'2.{x) and F'2,{x) are defined by 



and 



-3 



n\,n2=—oo 
+ 



) 

^niof + {n2h + [n26 + a;2]' + 

{riia + xif + {n2hf' [n2hf + 
(nio + xi)' + (^26 + X2)' [n2h + X2\ 



F2{x) 



J2 \{niaf + {n2bf 



-2 



+ 



"1,^*2 =— 00 



+ [(riio)' + (77,26 + ^2)' 



(nio + xi)' + (7126)' 



-2 



+ 



+ 



For i = 3 



We still need to calculate 



1 d d 1 

Y,—di(t>ndj(t>*n = lim— — C(s = o I 



dx^ dy"^ 
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For i=2 and j=l, we have 



47r2 ' 



where the function F2i{x) is defined by 



F2i{x) = £ [(nia)2 + M)'] [n26][nia] + 

- ^{n2b + ^2)^ + (riia)^ [7126 + X2] [nia 



+ 



+ 



(nia + xi)^ + (n26)^ [nia + Xi] [7126] + 
[(nio + xi)^ + (^26 + ^2)^ [riio + xi] [^26 + X2] 



^ ~ Y (^1^ + ^1)^ + (^2^ + 2:^2)^ ['T'la + xi] [n2b + X2] , 

ni,n2=—oo 

because the first three summands are odd in one index. 
For i=3 and j=l and for i=3 and j=2, we have 



1 1 



n 



n 



To obtain the components of the improved stress tensor we need to calculate 

Y —(t>ndU*n ^Y^ (^Un) K = I™ ( ^ ) C(s = ^ I X, y). 



For i — 1 
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where the function -Dii(x) is defined by 



^11 (^) = H {{niaf + {n2bf [riia 



ni,n2 =—c>o 



+ 



(?T,ia + a;i) + (r22&) [uia + xi] + 
(nia + Xi)^ + (?22& + 2:2)^ [nio + xi]^} 



(53) 



For i = 2 



(54) 



where the function .022(2;) can be obtained from Dii{x) with the change Xi ^ X2 and nia ^ 77.26. 



For 2 = 3 



1 



1 



167r2 



Fo(f). 



(55) 



Finally 



For i = 1 and j = 2 



1 9 9 1 

^ —^ndidjcpl = lini — — C(s = - | x,y). 



Fi2{x) 
47r2 



Substituting the results of eqs.(^, (gD|), (|D, 



(56) 



(57) 



and ( pop in eqs. ([T^-|TB|), we obtain: 



(roo(f)) = -^(Di(f)+D2(f)) + ^(Fi(f)+F2(f)) 



(Too(f))^ + (Too(f)) 



F ' 



(5^ 
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where B and F mean the boundary divergent and finite part respectively. In exphcit form: 



(^oo(^))i 



327r2 



+ 



[a - xif + {b- 



+ 



[a - xi) + {X2) } 



+ 

-2, 



167r2 



{[X2] ^ + [6 - X2] ^ + [xi] ^ + [a - xi] 



(59) 



327r2 



+ 



(ni,n2)7^(0,0) 



E 



+ 



(nia + xi) +(n26 + X2) 

(ni,n2)^(0,0), (-1,-1), (0,-1), (-1,0) 

+2 [{niaf + {n2h + X2fY^ {[n2h + X2f - [niaf) + 

(ni,n2)7^(0,0),(0,-l) 

+2 E 

(ni,n2)7^(0,0),(-l,0) 



2 2 " — ^ 2 2 

(nia + xi) + (n26) ([nia + xi] -[^26])). (60) 



The restriction in the first sum above accounts for the exclusion of the free space divergent term. 
We see that {Tqq{x))q diverges in all walls, i.e., xi = 0, a and X2 = 0, 6 and in all edges, (xi, X2) = 
(0, 0),(a, 0),(0, 6), (a, h) of the waveguide. The other components are: 



(r33(f)) = - (roo(f)) 

(Ti2(f)) = (T2,(a;)) : ^ 



87r2 



Fl2(^) 



(T23(f)) = (T32(f)) = (T3i(f)) = (Ti3(a;)) = (To,(f)) = (T,o(a;)) = 0. 



(61) 

(62) 
(63) 
(64) 
(65) 
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It can be verified by writing these quantities explicitly that all the non-zero components have 
boundary divergences in all walls and all edges also. As has been remarked previously by many 
authors the divergences that appear in some components of the vacuum expectation value of the 
stress tensor are related with the unphysical boundary conditions imposed on the field. We can 
understand why the renormalized stress tensor becomes infinite on the boundary. This is related 
with the uncertainty relation between the field and the canonical conjugate momentum associated 
with the field ITTI, ITBlira. 



For the improved stress tensor we have 



(eoo(x)) 



1 



487r2 



= - (Too(x)) ^Fo(x) 



(66) 



Explicitly: 



(eoo(f)) 



B 



967r2 



-2 



+ 



[a - xif + (6 - X2) 



-2 



+ 



(xif + {b- XiY 



[a - xif + {x2)'^ } 



(67) 



(0oo(^))i 



+ 



3277 

1 

24^ 



1 —2 1 

— ^ Uriiaf + {n2bf] - —— ^ \{n,a + Xi^ + {n^h + X2)' 



^{n2hY [riia + XiY + {n2l)f'\ +^{niaf' {niaf + {n2h + X2 



(68) 
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We have to exclude: in the first sum the term (ni,n2) = (0,0), in the second sum the terms 
(ni,n2) = (0,0), (-1, -1), (-1,0), (0, -1). 

We see that (6oo(a;)) has no wall divergences but only edge ones, as pointed in ^. Although 
we cannot associate a curvature length to the edges, they seem to have a similar behaviour, since 
the divergences associated with them still remain even for the conformally coupled scalar field. 
The other components are: 

(0n(.l) - -^^DM + ^l^.(f) - ^l^n(f) + ^F,(x) - ^F.(x) (69) 

(e33(x)) = - (eoo(x)) (71) 



(ei2(f)) = (621(5)) = -^^12(5) (72) 



(623(5)) = (032(5)) = (631(5)) = (9x3(5)) = (00.(5)) = (e,o(5)) = 0. (73) 

Again, by writing these quantities explicitly, it is easy to see that all non-zero components of the 
improved stress tensor are free of wall divergences but have edge divergences. 

We are now interested in comparing the local calculation of the energy density (Too (5)) with the 
more familiar global one. In this way, let us now calculate the global energy inside the waveguide 



by integrating the energy density (Too (5)) given by eg. (^8]) over the cavity, < xi < a, < X2 < &. 
Despite the fact that a closed form of the double sums in (Too (5)) are presently not known, it is 
possible to calculate its integrals over the spatial region of the waveguide. We shall devide this 
total energy by the area of the cross-section of the waveguide a x b, and it is usually refered to 
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also as energy density, although this comes from an integrated quantity per unit area and is not 
actually a density in the sense of a local quantity, this one legitimately represented by (Too(x)). 
We shall assume, for the global computation, that the field exists only inside the cavity. 
Clearly the following expression: 



dxi / dx2 {Too{x)) = {Too{x)) 

Jo J J cavity 



cavity 



{Tooix))^+ [ [ (roo(f))^ (74) 

J J cavity 



diverges because of the first term: {Tqq{x))^ is divergent on the walls and edges. So let us treat 
the second term, in which {Tqq{x))p is given by eq. (|60D and is finite. The integral of the first term 
of eq. (|60D gives: 



ab 
32^ 



Z{2\a,b), 



(75) 



where, in the notation of 



Z{2\a,b)= J2 [M^ + {n2b)' 

{ni,n2)7^{0,0) 

Again following the notation of [0, the integral of the second term of eq.(|60|) is given by: 



(76) 



1 r 



327r2 



pa no 

I dxi / (ix2 CF(2|a;i, X2) 
Jo 



327r2 



Cb(2|xi,X2) 



87r2 



ext 

00 roo roo rb pa poo 

/ + / / + / / 

a Jb Ja Jo Jo Jb 



dxdy 



(x^ + y'^Y ' 



(77) 



where 



Ci.(2|xi,X2)= E 

(rii,n2)^{0,0)(0 -1)(-1,0)(-1,-1) 



(nia + x{f + (77-26 + X2)' 



(78) 
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and the right-hand side is each one of the divergent edge terms integrated over the appropriate 
quadrant outside the cavity, i.e., away of the points where they diverge, and thus eq.(^) is also a 
finite contribution (see for further explanations). 

The integration of the third (and fourth) term of eq.(pD|) is not difficult: 
1 



^r, o I dxi dx2 
IbTT'^ Jo Jo 



E 



(ni,n2)^(0,0)(0,-l) 



(riia)^ + (^26 + X2)^ (^26 + X2)^ - (riia)^ 



-Tib r^'^-'X"^-^ 



+ 



+ [{riiaf + {n^b + X2y] [(«2& + 2:2)^ - (^la)^] ] • (79) 

n2=—oo,ni^0 



The first integral above was also calculated in [0: 

r / dxi [ dx2 7 — ; — 

16n^ Jo Jo „2^-i (^20 4- 



+ X2 



ext 



Cb(2|0,X2), 



where 



Cb(2|0,X2) = 

ext Jb 



1 /-o 

X2 ^ —00 



dx2 —+ I dX2 



{b - X2y 



(80) 
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is also finite. The other term gives 
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where 



00 

C(s) = ^n- 

n=l 
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is the usual Riemann zeta function, and use has been made of the integral K 







[1 + px\ V[y) 



Gathering all previous results we have that: 



where 



/ (roo(f))5 = --^ (a / CBmx2) + b[ C^(2|xi,0)) / CB(2|xi,a;2) (84) 

Jext lOTT^ V Jext Jext J SZTT^ J ext 

is finite, because it is the sum of each of the wall and edge divergent terms integrated outside the 
cavity, i.e., far from the spatial points where they diverge. Eq.(^3|) can be written as: 

^ f f (Too(f))^ = Ec{a, b) + ^ f {Too{x))^ , (85) 

do J J cavity OO Jext 

where Ec{a, b) is the global Casimir energy divided by the cross-section area ax 6 for the waveguide, 
in agreement with (in fact. Actor's definition of V^jf is twice the usual one). We can add the 
same infinite term 

(Too(f))B 



/ cavity 

to both sides of the equation above, obtaining: 

\ f f (roo(x)) = Ec{a,b) + - ( {Too{x))b, (86) 

0,0 J J cavity OO Jail space 

where the last integral above is an infinite constant independent of the cavity dimensions a, b. 
In global calculations one usually discards this infinite constant because it does not give rise to 
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forces. Discarding this infinite constant, one obtains from the expression above the total Casimir 
energy per unit area inside the waveguide with Dirichlet boundary conditions in all walls. It can 
be shown that the improved stress-tensor yields the same Casimir energy per unit area Ec{a,b): 

^ ^ ^ (eoo(x)) = Ec{a,h) + \ ( (87) 



ah J J cavity ah . / all space 

It is known that the sign of the global Casimir energy is dependent on the relative size of a 
and h. For example, for the square waveguide a = 6 a positive value for Ec{a, h) is found. Because 
this is a symmetric configuration, an equal total outward force appears acting on each of the four 
walls, which tends to make the cavity expand. 

An important lesson that we learn from eq. (P5|) is that the integral inside the waveguide of 
the finite part of (Too(s)) does not yield directly the total energy Ec{a,h), but this one plus the 
constant: 

C{a,h) = 1/ (Too(f))^- 



ext 

1 
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Dowker and Kennedy have evaluated the total energy of the conformally coupled scalar field 
in the interior of the waveguide for two special configurations. For the square a = h, they showed 
that it assumes a positive value. When h = 2a the energy decreases, assuming a negative value. 
Figures ([I|) and (|^) show the form of {Tqq{x))p and {Qoo{^))f square waveguide, assuming 

a = h = 1. They present a minimal value in the middle of the waveguide and assume only positive 
values, which produces a positive value. From the integral of this density one should subtract the 
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constant C(l, 1) in order to obtain the total Casimir energy per unit area of the square waveguide. 
As the value of b increases (for a = 1), these local quantities acquire negative values in some space 
points, making the total energy decrease. Figures (|^) and (|D show the local energy density in 
the case b = 2a. In this case the contribution of the negative part of the local energy dominates 
and since one still has to subtract C(l,2) from the integral of this energy density, one obtains a 
negative total energy per unit area. 

4 Local forces 

In this section we will calculate the local Casimir force density that acts on the walls of the 
waveguide. To do this we will use the relation between the local force density and the discontinuity 
of the stress tensor across the walls. Although we don't know the modes outside the waveguide 
(because the external mode problem for the waveguide is unsolved), we can introduce an external 
structure where the modes of the field are known 0, in such a way that the interior region is the 
interior of the waveguide. One way to do this is connecting two parallel infinite Dirichlet planes 
by two strips. In this configuration, we know the modes in all regions and the stress tensor can be 
calculated anywhere. Let us position two parallel infinite Dirichlet planes at Xi = and Xi = a 
and connect these planes by two strips, positioned at X2 = and X2 = b. The interior region of this 
configuration is just the waveguide. In the regions xi > a and Xi < there are no contributions 
from the stress tensor to forces that act in the two infinite planes {{Tu{x)) = 0). In the regions 
< Xi < a, X2 > b and < Xi < a, X2 < 0, the components of the stress tensor have a nonzero 
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contribution to the forces that act on the strips. In these regions the stress tensor has already 
been calculated in p. For completeness, we present the relevant component here, i.e., (T22 (a^))exr 



{T22 m 



ext 
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-3 



[na + Xi)^ + X 



+ 4x0 



1 + 



na + Xi 

X2 



(89) 



The equation above will serve to compute the local Casimir force that acts on the strip at X2 = 
(a similar one exists for the strip at X2 = b). We note that the edge divergences above at 
(xi, X2) = (0, 0), (a, 0) will not be canceled, when we come to calculate the local force, by those of 
the interior of the waveguide that appear in eg. (|62D . Nevertheless neither the equation above nor 
eq.(^) present wall divergences as X2 — > 0. Thus the local Casimir force at the strip at 0:2 = 
diverges only at the edges, but not on the strip. 

We note also that the components (T21 {x))^^^ and (T12 {x))^^^ vanish on the walls. To obtain 
the local forces, we use the local force density that acts on the point x and is given by fi (x) = 
—djTij (x) . Thus the local force per unit area on the boundary plane at Xi = is: 

^^""'^ - limJ(Tn (xi = -£))- (Tn (xi = e))] 



1 



327r2 



n-i ,n2=— 00 



'4:^{niaY + {n2b + X2Y [niaf + 



{niaY + {n2bf 



{uiaY + {n2b + X2f ), 



(90) 



in the positive Xi— direction, and an equal but opposite force acts in the plate at Xi = a. (We have 
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to exclude the term {rii = n2 = 0) in the third and fourth sums and the term {ni = 0, ^2 = — 1) in 
the fourth sum. The last two exclusions accounts for the renormalization of the edge divergences.) 
The force on the wall parallel to the plane at X2 = is given by 



A 



hm [(T22 {X2 = -e)) - (T22 ix2 = e))] 
00 

2 {{any^ - {an + xiy^} + 



1 



327r^ 



1 



327r2 



ni,n2=—oo 



+4[{nia + xif + {n2b)^] ^ [n2bf + 

+ ^{riiay + {n2by - {uia + Xi)"^ + {n2b)'^ ) 



1 



327r2 
1 



E 



-3 



{niay + {n2by [^26] - 4 {uia + XiY + (n2&)^ ^fe] + 



+ 



167r2 



,n2=— 00 
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(nia + xif + {n2by 



(91) 



ni=— 00 712=1 

An equal but opposite force acts on the wall at X2 = b. (The edge divergences of eq.(|89D do not 
cancel those of eq. (|52D , as we have stressed; nevertheless these were discarded when calculating 
the local force above, and thus they do not appear.) 

Let us analyse how the local forces calculated above depend on the relative sizes of the waveg- 
uide. Figure (|^) shows the dependence on X2 of the finite part of the local force that acts on the 
wall parallel to the plane Xi = 0. It assumes only negative values and thus it is a repulsive force, 
in agreement with global calculations. The modulus of the force has a minimum in the middle 



of the wall and two maxima near the edges. Figure (Bl) shows the depence on xi of the force on 
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the wall parallel to the plane X2 = 0. It is an attractive force but with only one maximum in the 
middle of the wall. Although the global computation for the square waveguide gives a repulsive 
force in all walls, our attractive result is due to the external structure. 

Figures (0) and show the forces that act on the walls at xi = and X2 = when b = 2a. 
The local force at xi = assumes only positive values which makes it an attractive force, as we 
expect by approaching the parallel plate configuration, but still highly non-uniform. The force at 
X2 = assumes only positive values and it is very small in comparison with the previous force. As 
b grows, this force vanishes and the force at xi = behaves like the uniform Casimir force in the 
parallel plate configuration as figure shows. 

5 Conclusions 

In this paper we obtained the canonical and the improved stress-energy tensors of a massless scalar 
field in the interior of an infinitely long waveguide. The result found is strongly position dependent 
as expected. Although the global Casimir effect is related to experiments where we measure 
the force between macroscopic surfaces, the local properties of the vacuum field fiuctuations can 
in principle be observed by measuring the energy level shift of an atom interacting with the 
electromagnetic field. In the case of the local problem, surface and edge divergences appear 
related with the uncertainty principle. In order to compute the local forces we introduced an 
external configuration for which it is possible to solve the eigenmode problem. We have shown 
that the particular external configuration that we chose was not able to eliminate the wall and 
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edge divergences of the interior of the waveguide. In order to ehminate them two possible ways 
are to take into account the real properties of the material, i.e., imperfect conductivity at high 
frequencies, or else make a quantum mechanical treatment of the boundary conditions, as was 
done by Ford and Svaiter An alternative method of calculation (using a modified version of 



the Green's function method) to find the renormalized stress-energy tensor associated with the 
scalar field defined in the interior of an infinitely long waveguide is under investigation by the 
authors. 

We have also shown that the integral inside the cavity of the local result gives the known 
values for the global calculations, although the integral of the finite part of (Too(x)) gives the total 
Casimir energy plus a constant dependent on the waveguide sizes C{a,b). 
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1 1 



Figure 1: Renormalized local energy density of the minimally coupled scalar field in the interior 
of the square waveguide. 



31 



1 1 



Figure 2: Renormalized local energy density of the conformally coupled scalar field in the interior 
of the square waveguide. 
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2 1 



Figure 3: Renormalized local energy density of the minimally coupled scalar field for the b = 2a 
waveguide. 
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2 1 



Figure 4: Renormalized local energy density of the conformally coupled scalar field for the b = 2a 
waveguide. 
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Figure 5: Renormalized local force density that acts on xi = wall for the square waveguide. 
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Figure 9: Renormalized uniform force density that acts on xi =0 wall when b >> a. 
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